1 



On 
On 



O 

NO 



Non-perturbative improvement of operators with Wilson fermions* 

C. Dawson a , G. Martinelli b , G.C. Rossi c , C.T. Sachrajda a , S. Sharpe d , M. Talevi c and M. Testa b 

a Dept. of Physics, Univ. of Southampton, Southampton S017 1BJ, UK 

b Dip. di Fisica, Univ. di Roma "La Sapienza", P.le A. Moro 2, 1-00185 Roma, Italy 

c Dip. di Fisica, Univ. di Roma "Tor Vergata", Via della Ricerca Scientifica 1, 1-00133 Roma, Italy 

d Physics Dept., Univ. of Washington, Seattle WA 98195, USA 

c Dept. of Physics & Astronomy, Univ. of Edinburgh, The King's Buildings, Edinburgh EH9 3JZ, UK 
We outline two methods of constructing improved composite operators using Wilson fermions. 
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1. Overview 

A major source of errors in present lattice cal- 
culations is the use of a finite lattice spacing, a. 
This is particularly true for the calculations of 
properties of hadrons containing c and b quarks, 
which contain discretization errors proportional 
to am c and amb, respectively. Since lattice cal- 
culations have the potential to make significant 
contributions to the study of D and B physics, it 
is important to control such errors. 

One approach to solving this problem is the 
Symanzik improvement program. This allows one 
to remove discretization errors order by order in 
a. This program has been pursued by the AL- 
PHA collaboration, who have used chiral Ward 
Identities to determine, non-perturbatively, the 
on-shcll improved action, the improved vector 
current including am terms, and the improved 
axial current, scalar and pseudoscalar densities 
in the chiral limit Q-q]. Ward Identities for 
non-degenerate quarks allow the determination of 
some of the 0(am) improvement terms, at least 
in the quenched approximation Q. These meth- 
ods do not, however, allow one to determine the 
improved pseudoscalar density, axial current or 
tensor density away from the chiral limit. 

To do so we have developed two methods which 
remove all errors of 0(a), including those of 
0(am), from the matrix elements of bilinears. 



These involve (1) matching correlators of bilinears 
to their continuum form at short Euclidean dis- 
tances, and (2) matching quark and gluon corre- 
lators to their continuum form at large Euclidean 
momenta. Both rely on the restoration of chi- 
ral symmetry at short distances, and both work 
for quenched and full QCD. The first method re- 
quires only on-shell improvement, i.e. improve- 
ment of physical quantities, while the second re- 
quires improvement of off-shell correlators at an 
intermediate stage. 2 Details of the first method 
are given in Ref . H . Details of the second method 
will be forthcoming, along with results of a nu- 
merical pilot study |6| . Here we provide a sketch 
focusing on some important issues. 



2. Improvement Program 

We need first to improve the action Jl[]. 
continuum notation, the Wilson action, 



S 



w 



2y; 



Tr(F^F^) + qty + m )q 



In 



(1) 



*Talk presented by S. Sharpe. 



has 0(a) errors due to the derivative. Improve- 
ment requires the addition of all dimension 5 op- 
erators allowed by the symmetries, 

+ b g amTr(F f _ l „F filJ )/(2gl) -b m am 2 qq 

2 Throughout we use "improved" to denote quantities in 
which all errors proportional to a (multiplied by any power 
of In a) have been removed. Errors of 0(a 2 ) remain. 
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+ c[ amq(ft + m Q )q + c' 2 aq(ft + m ) 2 q,(2) 

with appropriately chosen, go-dependent, coeffi- 
cients Q. Here m = too — mdgo) is proportional 
to the physical mass, with m c the critical mass. 
The role of the various terms is as follows: 

1. The csw term improves dimensionless physi- 
cal quantities. Its determination is clearly essen- 
tial, and has been carried out in Ref. 

2. The coefficients b g and b m remove 0(am) 
terms from the relation between bare and renor- 
malized couplings and masses, respectively. For 
example, b g introduces a mass-dependence in the 
effective gauge coupling, g 2 s = g^iX — b g am), in 
such a way that the lattice spacing remains fixed 
as one varies to at fixed go- We will need to de- 
termine b g in our first method ||. 

3. The coefficients and c' 2 are only needed for 
off-shell improvement. One way to see this is to 
note that they can be removed by a change of 
quark variables (and an appropriate shift in go). 
This affects external sources, but not the spec- 
trum. In fact, one can ignore these terms also 
when doing off-shell improvement, because they 
can be absorbed by a suitable change in the im- 
provement coefficients for quark fields Q . 

The next step is to improve the operators them- 
selves. On-shell improvement has been discussed 
in Ref. 0; off-shell improvement requires addi- 
tional terms. For example, the improved form of 
the bare pseudoscalar density P = qjaq is (for 
degenerate quarks) 



P{x) = Z P {g 2 , [ ia)(\ + bpam) (P(x) (3) 

+ac' P q 7 5 (^ + to ) + (-^ + to )75 q(x)^j . 

In addition to the two on-shell improvement co- 
efficients, Zp and bp, there is new coefficient, c' p . 
This multiplies an operator which vanishes by the 
equations of motion, and so does not contribute 
to on-shell matrix elements. The pattern is the 
same for other bilinears — each has a single addi- 
tional off-shell term H. 

3. Gauge-Invariant Method 

This method involves only on-shell quantities, 
and so we do not need consider the off-shell coef- 



ficients such as c'p. We discuss the example of the 
pseudoscalar density — details for the other bilin- 
ears can be found in ||. 

To determine Zp and bp we require that the 
Euclidean two-point function of the improved lat- 
tice operator (evaluated using the on-shell im- 
proved action), 



G P (x) = (P(x)P(O)) 



= Zp 2 {l + b P am) 2 {P(x)P(0)) , 



(4) 
(5) 



agrees with the continuum result up to 0(a 2 ). 
The latter can determined, at short distances, us- 
ing the OPE 

— 2a s jp hi(x[i) + . . . 



Gcont 
p 



(X) 



2ir 4 x 6 



(G) 



x [1 + O(toV) + 0(toA^ cd x 4 , A 



l QCD J 



')] 



The term on the first line is the coefficient func- 
tion of the unit operator, with \x the renormaliza- 
tion point and "fp the one-loop anomalous dimen- 
sion. The correction terms on the second line are, 
respectively, the perturbative contribution due to 
the quark mass, and the non-perturbative terms 
due to the operators ipi/) and F 2 respectively. The 
crucial point is that chiral symmetry is restored 
at short distances (i.e. there is no to dependence 
if power corrections can be ignored). 

Equating G p and G c p nt in the chiral limit yields 
Zp in the chosen renormalization scheme. De- 
manding that Gp contain no terms linear in the 
quark mass determines bp, i.e. 



(1 + bpam) 2 = Gp(x; to = 0)/Gp(x; to) , 



(7) 



where Gp is the bare lattice two-point function. 
We emphasize that the determination of bp is 
non-perturbative, and independent of that of Zp. 

There is one subtlety in the determination of 
bp. The condition (Q) requires that the physical 
distance x not depend on to. This in turn requires 
the determination of b g , as discussed above. Note 
that this is only an issue for full QCD b g = in 
the quenched approximation. To determine b g 
non-perturbatively one must hold fixed a phys- 
ical quantity which itself has no dependence on 
to. The choice we suggest is the force between a 
heavy q — q pair at short distances. See Ref. || 
for details. 
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The determination of Zp requires that there ex- 
ist a window at short distances: a <C |x| <§C AqJ, d . 
This ensures the smallness of both higher or- 
der discretization errors of 0(a 2 /x 2 ) and non- 
perturbative contributions. In addition, to de- 
termine bp, one must work in a region where 
m 2 x 2 <C ma and toAq CD 2; 4 <C ma, so that 
the dominant mass dependence is due to the mo 
terms. These conditions can be satisfied for suf- 
ficiently small a, and for a suitable choice of 
m. Based on experience with non-perturbative 
renormalization on quark states, we expect that 
a = 0.05 — 0.1 fm will be small enough to ap- 
ply the method, but numerical tests are clearly 
required to check this. 

A potential practical problem is the need to 
have both short and long distances on the same 
lattice. Short distances are needed for the deter- 
mination of the improvement coefficients, while 
long distances are required in order to calculate 
matrix elements using the improved operators. 
This problem is not particular to our method — 
to calculate renormalization constants such as Zp 
requires making contact with perturbation the- 
ory and thus working at short distances. In the 
Schrodinger functional method, this problem is 
overcome by matching between lattices of differ- 
ent spacings. We wish to stress that a similar 
approach is possible here. In essence one cal- 
culates the improvement coefficients with a very 
small lattice spacing, and then matches these co- 
efficients onto the lattices with larger spacing on 
which the matrix elements are calculated. The 
details are sketched in 

4. Gauge Non-invariant Method 

Our second idea is an extension to 0(a) of 
the non-perturbative renormalization program of 
Ref. 0| . We require that quark and gluon cor- 
relators agree with their renormalized continuum 
counterparts at large Euclidean momenta. Two 
major complications arise in this extension. First, 
since we are improving off-shell quantities, we 
need to include the additional improvement co- 
efficients such as c! P . Second, since we must fix 
the gauge, the improvement terms are no longer 
constrained by the gauge symmetry, but rather 



by the lattice BRST symmetry. This allows addi- 
tional improvement terms which are gauge non- 
invariant (or non-covariant). The only such term 
which appears at 0(a), however, is a (fq term in 
the improved quark field: 

q = Z q (l + b q am) [l + ac q (ft + mo) + acMGi<f\ q ■ 

Gauge invariant bilinears require non-invariant 
improvement terms only at 0(a 2 ). 

We can determine all the on-shell and off-shell 
improvement coefficients by requiring that chi- 
rality violating form factors vanish at large Eu- 
clidean momenta. For example, the amputated 
vertex of the improved pseudoscalar, P, should, 
for large momenta, be proportional to 75 with no 
a (lfi ~ 1/2)15 t erm - It turns out that the off-shell 
terms play an essential role in this procedure: one 
must first determine the off-shell coefficients in 
order to correctly determine all the on-shell co- 
efficients. A numerical test of this approach is 
underway. A similar method has been suggested 
in Ref. ||, although without the inclusion the 
gauge non-covariant cngi term. 

5. Conclusions 

We have proposed two types of improvement 
condition, both non-perturbative, and sketched 
their use for bilinears. They provide an alterna- 
tive to Ward Identities for determining improve- 
ment coefficients in the chiral limit, and they have 
the advantage of working also away from the chi- 
ral limit. This is particularly important for ap- 
plications involving heavy quarks. The meth- 
ods should be straightforward to generalize to 
more complicated composite operators. We are 
presently studying their numerical efficacy. 
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